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Abstract. We describe the interaction pattern in the x-y plane for a family of soliton solutions 
of the Kadomtsev-Petviashvili (KP) equation, 

(— 4m< + Ujcxx + 6uu x ) x + 3u yy = 0. 
Those solutions also satisfy the finite Toda lattice hierarchy. We determine completely 
their asymptotic patterns for y — > ±°°, and we show that all the solutions (except the one- 
soliton solution) are of resonant type, consisting of arbitrary numbers of line solitons in both 
aymptotics; that is, arbitrary N- incoming solitons for y — » — °° interact to form arbitrary N + 
outgoing solitons for y — > °°. We also discuss the interaction process of those solitons, and 
show that the resonant interaction creates a web-like structure having (JV_ — — 1) holes. 
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1. Introduction 



In this paper we study a family of solutions of the Kadomtsev-Petviashvili (KP) equation 

d_ 

dx 

which can be written in the bilinear form 1 7 1 



dt 3x 3 dx 



d 2 u 



(1.1) 



[ - AD X D, + D A X + 3D] ) x • x = . 



(1.2) 

Here D x ,D y andD, are the Hirota derivatives, e.g., Dff ■ g = (d x — d x i) m f(x,y,t)g(x',y,t)\ x=x i 
etc., and u is obtained from the tau-function x(x,y,t) as 



d 2 

u(x,y,t) =2^- 1 \ogx(x,y,t). 
ax L 



(1.3) 



It is well-known that some solutions of the KP equation can be obtained by the 
Wronskian form T = 1m (see Appendix and also |4|), with 

/i0O ••• f M 0O 



T M = Wr(/i,...,/ M ) 

' fi0M-l ■■■ f M 0M-l 

§ To whom correspondence should be addressed (kodama@math.ohio-state.edu) 



(1.4) 
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. ,M} is a linearly independent set of M 



for 1 <i <M. 



(1.5) 



where f0n = d"f/dx", and {f{x,y,t) \ i = 1 
solutions of the equations, 

df = ^f df L= #fi 
dy dx 2 ' dt 8x 3 

For example, the two-soliton solution of the KP equation is obtained by the set {f\,f2}, with 
/. = e e 2i _ 1+e e 2i) i = i,2, (1.6) 

where the phases ; are given by linear functions of (x,y,t), 

Bj(x,y,t) = -kjX+kjy-k 3 jt + Q j, j=l,...,4, (1.7) 
with k\ < Ic2 < < k\. This ordering is sufficient for the solution u to be nonsingular. (The 
ordering k\ ^ ki < kj ^ k/\ is needed for the positivity of T2.) Note, for example, that if 
k\ < ki < k2 < k\, 12 takes zero and the solution blows up at some points in (x,y,t). The 
formula il.6\ can be extended to the M-soliton solution with {f\ , • • • 1 8 1. 

On the other hand, it is also known that the solutions of the finite Toda lattice hierarchy 
are obtained by the set of tau-functions {xm \ M = I,... ,N} with the choice of /-functions, 

N 

fi=Le 9i =:f, 

i=i (1.8) 

fi=f0i-l, l<i<M<N, 

where the phases 0,, 1 <i <N are given in the form il.7\ (see for example 1 13 1). This implies 
that each tau-function Xm gives a solution of the KP equation. If the /-functions are chosen 
according to i ll ,8t . the tau-functions are then given by the Hankel determinants 

/0O ••• /0M-1 

: ■.. : , for I <M<N. (1.9) 

f0M - 1 • • • f02M - 2 

Note here that = C exp(9i H h Qn), with C = constant, yielding the trivial solution, and 

1m and 1n-m produce the same solution with the symmetry (x,y,t) — > (— x, —y, —t), due to the 
duality of the determinants (i.e., the duality of the Grassmannians Gi(M,N) and Gr{N — M,N); 
see also Lemma lzTY The finite Toda lattice hierarchy is defined in the Lax form 1 3 1 
dL 



dt„ 



l,...,N-l, 



(1.10) 



where the Lax pairs (L,B n ) are given by 



B„ 



1 bi fli 












a\ b2 


a\ 








«2 















bN-i 


OiV-l 




V 





on -1 


b N 


) 


|((^)>o- 


-(L" 


<o) , 







and where C>o (C<o) denotes the strictly upper (lower) triangular part of a matrix C. Here the 
flow parameters ?,'s are chosen as t \ = x, t2=y and 1 3 = t for the KP equation. The functions 
[ai,bj) are expressed by 



On a family of solutions of the KP equation satisfying the Toda lattice hierarchy 



3 



where Xq 



1 . Then the tau-functions x„ satisfy the bilinear equations 
1 



-D 2 l -T 



T11T17 



(Xn,; 



: Xn+lXn-1 ■ 



(1.12) 



which are just the Jacobi formulae for the determinants D 

1 
1 



-1. i-e. 



D 



D 



D 



n 

n+l 



D 



n+l 
n 



D 



n, n 
n, n - 



D . 



Here D 



hi 
k,l 



(1.13) 

denotes the determinant obtained by deleting the z'-th and j-th rows and the 



Ic-th and the l-th column in D \ 5\. 

Remark 1.1. According to the Sato theory (see for example 1 14 1), these bilinear equations for 
the KP equation and the Toda lattice hierarchy are the Pliicker relations with proper definitions 
of the Pliicker coordinates ly labeled by Young diagrams Y = {l\ , £2), with l\ < £2 giving the 
numbers of boxes in Y, 

T (0,0)X(2,2) - t(0,l) x (M) + T (0,2)X(1,1) =0. (1.14) 

For the KP equation, those Pliicker coordinates are related to the derivatives of the tau-function 



T (0,0) 


= Xm, 




X(0,1) 


= d x X M , 




X(0,2) 


= Wh 


"3y)XM, 


T (l,l) 




" 3y)XM 


T (l,2) 




-3;)Xm, 


. T (2.2) 


=m 





Then the Hirota bilinear equation (II. 2> is equivalent to the Pliicker relation J 1 . 1 4b . For 
the Toda lattice equation, the Jacobi formula Jl . 1 3i can be considered as (I1.14> with the 



identification x 



(0,0) 



D etc. 



as x - 
as x 



(1.15) 



We should also recall that the solutions of the Toda lattice equation show the sorting 
property of the Lax matrix L \ 12 1; that is, 

diag(X,i,...,A,jv) 
diag(Xjv,...,A,i) 

where X\ > X2 > ■ ■ ■ > Xn are the eigenvalues of L. These eigenvalues are related to the 
parameters in (II. 7> as A,; = —hi (see below). 

In this paper, we are concerned with the behavior of the KP solutions Jl .31 whose tau- 
functions are given by II. 9> . We describe the patterns of the solutions in the x-y plane where 
each soliton solution of the KP equation is asymptotically expressed as a line, namely, 

x = c±y + £± for y — > ±°o 

with appropriate constants c± and f;± for a fixed t. In particular, we found that all the 
solutions (except the one-soliton solution) are "resonant" solitons in the sense that these 
solutions are different from ordinary multi-soliton solutions. The difference appears in the 
process of interaction, which results, for example, in a different number of solitons (or lines) 
asymptotically as y — > °o or y — > — 00, 
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In our main result (Theorem 12. 51 we show that for the solution with the tau-function 
given by II. 9> with 11.81 . the number of solitons in asymptotic stages as y — > ±°°, denoted by 
N + and A^_, is given by 

N+=M, N-=N-M. 

Thus, the total number of exponential terms in the function / in J 1 .81 gives the total number 
of solitons present in both asymptotic limits, i.e., N = + N + , and the number of outgoing 
solitons N+ is given by the size of the Hankel determinant II. 9> . We call these solutions 
"(Af-,Af+)-solitons". In particular, if N — 2N+ = 2N-, the solution describes an A^-soliton 
having the same set of line solitons in each asymptotics for y — > ±°°. However, these multi- 
soliton solutions also differ from the ordinary multi-soliton solutions of the KP equation. The 
ordinary n-soliton solution of the KP equation is described by n intersecting line solitons 
with a phase shift at each interaction point. If we ignore the phase shifts, these n lines form 
(n — 1) (n — 2)/2 bounded regions in the generic situation. However, the number of bounded 
regions for the (resonant) Af + -soliton solution with II. 9> is found to be (A^ + — 1 ) 2 ; for example, 
even in the case of a two-soliton solution there is one bounded region as a result of the resonant 
interaction. In general, we show in Proposition 13. II that for the case of a (A r _,A r +)-soliton 
solution, the number of bounded regions (holes) in the graph is given by (N- — i){N + — 1), 
except at finite values of t in the temporal evolution. 

These resonant Af+-soliton solutions are similar to some of the solitons of the coupled KP 
(cKP) hierarchy recently studied in Ref. |9|, where such solutions were called "spider-web- 
like" solutions. The analysis of finding web structure that we describe in the present study 
may also be applied to the case of the cKP hierarchy. 



2. Asymptotic analysis of the solutions 

Before we discuss the general case for the tau-function ( II .91 with II. 8> . we present some 
simple cases corresponding to a (l,l)-soliton and a (2,l)-soliton solution; the latter turns out 
to be the resonant case of an ordinary 2-soliton solution of the KP equation. 

As explained in the Appendix, we first note that the (N-, 1) -soli ton can be described as 
the solution of the Burgers equation IA.3> . 

dw\ _ dw\ d 2 w\ . , d 

-5— +2wi^— = -z-j-, with wi = - — logTi. 

ay ox ox ox 

An explicit solution of this equation is a shock, which corresponds to the case ofA r = 2 = l + l, 
i.e., Ti = e 6 ' + <? e 2. The solution w\ is then given by 

wi = £(*i+* 2 ) + 3(*i-*2)tanh£(8i -6 2 ) 

k\ as x — > 00 j 
, (for ki < k 2 ) 

k2 as n —00, v ' 

which leads to the one-soliton solution of the KP equation, 

d 2 

u = 2-^logx l = ±(k l -k 2 ) 2 sech 2 ±(Q 1 -Q 2 )- (2.1) 

In the x-y plane, this solution describes a plane wave u = <b(k x x + k y y — cof ) having the 
wavenumber vector k = (k x ,k y ) and the frequency CO, 

k= (-h+k 2 , k\-k\) =: ki, 2 , co = k\ -k\ =: CO1.2. 

Here (k, co) satisfies the dispersion relation, 4co£ v + k\ + 3k 2 = 0. We refer to the one-soliton 
solution 12. II as a line soliton, which can be expressed by a (contour) line, 0i = 2 , in the x-y 
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Figure 1. The confluence of two shocks of the Burgers equation with {k\ , k% , £3) = (— ^ , j , 1), 
which also represents a resonant soliton solution of the KP equation. Here and in the following, 
unless indicated otherwise, the horizontal and vertical axes are respectively x and y, and the 
graph shows contour lines of the function u(x,y,t) = — 2d x wi (x,y,t). 



plane. In this paper, since we discuss the pattern of soliton solutions in the x-y plane, we refer 
to c = dx/dy as the velocity of the line soliton in the x direction; that is, c = indicates the 
direction of the positive y-axis. 

Now we consider the case of a (2,l)-soliton, whose tau-function is given by 

This situation is explained in Ref. 1151 . and the solution describes the confluence of two 
shocks. Taking k\ < k 2 < k^ without loss of generality, for y — > — °° the two shocks (which 
correspond to line solitons for u) have velocities c\ .2 = k\+k 2 an d c 2,3 = k 2 +£3, and the 
single shock for y — > °° has velocity C13 =k\ + k^. This case is illustrated figure [J with 
(ki,k2,k-yj = (— A simple analysis (see below for more details) shows that the 
function w\ = — d x \ogi\ takes the following asymptotic values: w\ ~ k\ = —\ for large x, 
w\ ~ kj, = 1 for large —x, and in the middle region for large —y, W\ ~ kj_ = i. 

This Y-shape interaction represents a resonance of three line solitons. The resonance 
conditions for three solitons with the wavenumber vectors {k,-,-| 1 < i < j < 3} and the 
frequencies {co,-,,- 1 1 < i < j < 3} are given by 

ki,2 + k 2 .3 = k K 3, and C0i, 2 + CO2.3 = ©13, (2.2) 

which are trivially satisfied by those line solitons. Here we point out that this solution is also 
the resonant case of the ordinary 2-soliton solution of the KP equation. As we mentioned 
earlier, the ordinary 2-soliton solution is given by the M = 2 tau-function dl.4> with II. 6> . 
The explicit form of the T2-f unc tion is 

T 2 = (*i - fc 3 )e ei+63 + (Jfci - k 4 )e^ +e * + {k 2 - h)e^ + (k 2 - k A )e^- +Q \ 

where, as before, 0,- = — k\X + kfy — k?t + Q®. Note that if k 2 = kj, the T2-function can be 
written as 

T2 = e e 1+ e 2 +e 4 ^ _ h) Ag -e 4 + ^ _ U)e -*i + (fe -k A )e^] , 

where A = exp(03 — 83) = constant. Since the exponential factor e e i+ e 2+ e 4 gives zero 
contribution to the solution u = 28^1ogX2, the ^-function is equivalent to the case of a (2,1)- 
soliton solution (except the signs of the phases, and more precisely it is a (l,2)-soliton); that 
is, the resonant solution with confluence of solitons. Note also that the condition ko = £3 
is nothing else but the resonant condition in Ref. II II . and it describes the limiting case of 
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an infinite phase shift in the ordinary 2-soliton solution, where the phase shift between the 
solitons as y — > ±°° is given by 

g= (h-h)(k 2 -k A ) 
{k2-kj,){k\ -fa) ' 

The resonant process for the (N- , 1 )-soliton solutions of the KP equation can be expressed as 
a generalization of the confluence of shocks discussed earlier. This case has been discussed 
in Ref. fT0l . 

We now discuss the general case of (N-,N+) -soli ton solutions. In order to describe the 
asymptotic pattern of the solution associated with the tau-function dl.9> . we start with the 
following: 

Lemma 2.1. Let f be given by 



N 

f=L 



with 9, = -kiX + kfy + Q?. 



Then for N = N+ + N- and 1 < N + < N — 1, the tau-function defined by the Hankel 
determinant M.9\ has the form 



= I A(h,...,i N+ ) exp £ di , 

i<h<-<i N+ <N \j=l 1 



(2.3) 



where A(i\ , . . . , iff) is the square of the van der Monde determinant, 
A(ii,...,i N )= n (kij-ki,) 2 . 

l<j<l<N + 



Proof. Apply the Binet-Cauchy theorem 1 5 1 for 



det 



e 6 ' 
he^ 



e Q 2 
fee 02 



k N e BN 



\k 



N. 



-1 



e" z 



( 1 fa 
1 k 2 



V 1 k N 



K N > 

N 



□ 



One should note from d2.3i that the -function contains all possible combinations of 
Af + phases from the set {8 ; - 1 j = 1, . . . , A^}, unlike the case of ordinary multi-soliton solutions 
of the KP equation. For example, the T2 -function for the 2-soliton solution with dl.6> include 
only four terms, and is missing the combinations 61+62 and 63 + 64. This makes a crucial 
difference on the interaction patterns of soliton solutions, as explained in this paper. In 
particular, we will see that the (A^_ ,N+) -solitons are all of resonant type in the sense that local 
structure of each interaction point in those solitons consists of either (2, 1)- or (l,2)-solitons. 

Remark 2.2. The % + -function given by fll .9I > is positive definite, and therefore the solution 
u has no singularity. In general, the Wronskian d 1 .4-1 takes zeros at some points in the 
flow parameters. The set of those points is called Painleve divisor, whose geometry has an 
interesting structure related to the Birkhoff stratification of the Grassmannian 1 1 1. Also, if 
one includes some exponential terms with negative coefficients in il Mi . the T^f-functions 
vanish on a set of points in (t\ , £2 , • • • , fjv- 1 ) ■ Then the set of those points can be described as 
intersections of the Bruhat cells of the flag manifold (see for example [2 1). 
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Figure 2. The functions T\j(c) = kj{kj-c) for (kxM-MM^h) = (-2,-l,-^,|,|). The 
levels of intersection 0, 1, 2, 3 are respectively denoted by circles, diamonds, triangles and 
stars. 

Let us now define a local coordinate frame (%,y) in order to study the asymptotics for 
large \y\ with 

X = cy + £ ) . 

Then the phase functions 9, in / of dl .81 become 

9i = -* i § + Ti i (e)y + e?, for i=l,...,N, 

with 

r\i(c) :=ki(ki-c). 

Without loss of generality, we assume the ordering for the parameters | i = 1 , . . . ,N}, 
k\ < lc2 < ■ ■ ■ < km ■ 

Then one can easily show that the lines r| = T];(c) are in general position; that is, each line 
T| = T[i(c) intersects with all other lines at N — I distinct points in the c-r| plane; in other 
words, only two lines meets at each intersection point. Figure |2] shows a specific example, 
corresponding to the values (^1,^2,^3,^4,^5) = (—2,-1, —A, |, |). 

Now the purpose is to find the dominant exponential terms in the Zn + -function (12.31 
for y — > ±°o as a function of the velocity c. First note that if only one exponential is 
dominant, then w\ — — 9^1ogTjv + is just a constant, and therefore the solution u — — 2d x w\ 
is zero. Then, nontrivial contributions to u arise when one can find two exponential terms 
which dominate over the others. Note that because the intersections of the r|,'s are always 
pairwise, three or more terms cannot make a dominant balance for large \y\. In the case of 
(N-, l)-soliton solutions, it is easy to see that at each c the dominant exponential term for 
y — > 00 is provided by only r|i and/or and therefore there is only one shock (W + = 1) 
moving with velocity c\ t N = k\ + k^ corresponding to the intersection point of rj 1 and r|^ 
(see figure 0. On the other hand, as y — + — °°, each term r| ; can become dominant for 
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some c, and at each intersection point r| ; = the two exponential terms corresponding 
to r\j and T| , + j give a dominant balance; therefore there are N- = N —I shocks moving with 
velocities cj = kj + kj + \, for j= l,...,N—l (see again figure[2}. 

In the general case, ^ 1, the Zn + -function in d2.3l > involves exponential terms having 
combinations of phases, and two exponential terms that make a dominant balance can be 
found as follows: Let us first define the level of intersection of T[i(c), 

Definition 2.3. Let r\i(c) and T\j(c) intersect at the value c = c/j = kj + kj, i.e., r|,(c,j) — 
T\j(cij). The level of intersection, denoted by (T/j, is defined as the number of other T|/ 's that 
at c = Cjj are larger than T\i(cjj) = T\j(cu). That is, 

C>U ■= Ml/ I ni(ci,j) > T\i(cij) = T\j(ci,j)}\ ■ 
We also define I(n) as the set of pairs having the level Ou = n, namely 

I(n) := {(r|;,r|;) | Lj = n, for i < j}. 

The level of intersection can take the range < O/.y <N — 2. Then one can show: 
Lemma 2.4. The set I(n) is given by 

I(n) = {(r|;,r|jv_„ +! -_i) I i = l,...,n + l}. 

Proof. From the assumption k\ < k2 < ■ ■ ■ < k^, we have the following inequality at 
c = Cjj (i.e. r\i = r\j) for i < j, 

r| /+ i , . . . ,riy_i < rii = i\j < tii, . . . ri ;+ i, . . . ,r\ N . 

Then taking j = N — n — 1 leads to the assertion of the Lemma. □ 

Note here that the total number of pairs (r|, ,r| ; ) is 



We illustrate these definitions in figure 13 where the sets I{n) for the level of intersection 
n = 0,1,2,3, which are respectively marked by circles, diamonds, triangles and stars, are 
given by 



Here recall that + = A^. These formulae indicate that, for each intersecting pair 
(T|;,r| ; ) with the level Af- — 1 (A^+ — 1), there are — 1 terms T|/'s which are smaller 
(larger) than r|, = T[j. Then the sum of those — 1 terms with either ri, or r| 7 provides 
two dominant exponents in the Xn + -function for y — > — °° (y — > °°) (see more detail in the 
proof of Theorem l2.5> . Note also that \I(N± — 1) | = . Now we can state our main theorem: 

Theorem 2.5. Let wi be a function defined by 




( 7(0) = {(Th,Tl5)}, 

I 7(1) = {(r)i,Ti4), (Tl2,r) 5 )}, 

I 7(2)={(T|i,ri3), (r| 2 ,Ti4), (r) 3 ,Ti 5 )}, 

{ 7(3) = {(T| 1 ,Tl2), (T| 2 ,ri3), (Tl3,T|4), (r) 4 ,Tl 5 )}. 



For the case of (A r _,A r +)-solitons, the following formulae are useful: 
f I(N--l) = {(m,r\ N++i )\i=l,...,N-}, 
\ I(N + -l) = {(m,T\ N -+i)\i=l,...,N + }. 



d 

W\ = -^logX7v + , 

with Xn + given by H2.31 . Then w\ has the following asymptotics fory — > ±> 



-oo 
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(i) For y — > — °° and x = c, \N + +iy + ^ far i — 1 , . . . ,N- , 

N++i 

K,(-,-):= L kj as $->-oo, 

iV + +i-l 

*(+,-):= I */ as 

(ii) For y — > °° and x = ci^+iy + £ /or / = 1 , . . . ,N+ , 

i-l 7V+-H-1 
«*(-.+):= I as 

i AT+-1 

*»(+>+) : = E I fev-y+i as £->oo. 

.7=1 .7=1 

where Cij = k[ + fry. 

Proof. First note that at the point T); = T[n + +U i- e -, { r [i, r [N + +i) S /(AL — 1), from Lemma l24l 
we have the inequality, 

1\i+l,1\i+2t ■ ■ ■ ,f\i+N+-l < T li = T l7V++i • 



N+-1 

This implies that, for c = kj + k^ ++ i, the following two exponential terms in the -function 
in Lemma im 

/JV + +i-i \ /iV++i 
exp(^ I 8,-J , exp^JE^fy 

provide the dominant terms for y — > — oo. Note that the condition T|,- = T|jv + +i leads to 
c = Cj^v++j = ki + kfj + +i- Thus the function w\ can be approximated by the following form 
along x = Ci !N++i y + 1, for y -> -«>: 

w, ^ log (A,(+, -) e -*iC + .-)5 + A,-(-, 

*,■(+, -)A,-(+, -)e-W+--K + £,-(-, -)A,-(- -y^j^ 
A / (+,-)e-^(+-^+A i (-,-)e-^(-H5 

£)(+, -)A,-(+, -) g fa++'-*^ + g f (-, -)A,-(-, -) 

A/(+,-)e (fciV ++ i ^ + A i -(-,-) 



where 



A,-(+,-)=A(j,...,iV + + z-l) expf £ 9° J 

A ( -(-, -) = A(i + 1, • • ■ ,N+ + i) exp ( 0° 

Now, from ki < kn + +i it is obvious that w\ has the desired asymptotics as \ — » ±°° fory — > — oo. 
Similarly, for the case of (Tjf, T[jv_+i) € — 1) we have the inequality 

Tl/='n2V_+i < Tli,Tl2,...,T|i_i,ri W _ + f + i ...,r\ N . 



N+-1 
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1.W+1 



N +1 W 



~Y °2,W+2 



C V 



w ,w 



"w ,w 



I J + "' + ^» I 

N +1 W 



'N -1.W-1 



/( w + , " + Vl 



2,W +2 



Figure 3. Asymptotic behavior of the function ivi with k\ <k% <■■■<§<•■■ < /cjv-i < k^. 
As y — » — °° there are /V_ jumps, moving with velocities Cj.ai + + ; (j = 1 , . . . , /V_ ) . As y — > oo 
there are /V + jumps, moving with velocities c;jv_ + ( (i = 1, . . . ,JV+). 



Then the dominant terms in the T/y + -function on x = q.al+O' + £ for y — * °° are given by the 
exponential terms 



exp( I Qj +\ l Qn-j+x ) , 

\j=i .7=1 



ex P ( L 8/ + L 6jv— y-t-l 



Then, following the previous argument, we obtain the desired asymptotics as \ — > ±°° for 

y — > oo. 

For other values of c, that is for c 7^ c,,/v + +; and c 7^ c, just one exponential term 
becomes dominant, and thus w\ approaches a constant as |y| — > 00. This completes the proof. 

□ 

Theorem 12.51 can be summarized in figure [3] As y — > — 00, the function wi has iV- 
jumps, moving with velocities cj^+j for j = l,...,iV_; as y — * oo, vvi has iV+ jumps, 
moving with velocities c, ^_ + , for i = 1, . . . ,iV+. Each jump represents a line soliton of the 
M-solution, and therefore the whole solution represents an (A^_,A^+)-soliton. Each velocity of 
the asymptotic line solitons in the (A^-,A^+)-soliton is determined from the c-r| graph of the 
levels of intersections (see figure|2j. For example, in the case of (l,4)-soliton in figure|2] one 
incoming soliton has velocity 01,4+1 = c\ 5, corresponding to the set 1(0), and four outgoing 
solitons have the velocities c^i+z for i = 1 , ... ,4, corresponding to 1(3). Note that, given a set 
of N phases (as determined by the parameters for i= 1 , . . . , AO, the same graph can be used 
for any (Af_,Af + )-solitonwithiV_+iV + =N. In particular, if N = 2M, we haveA^ + = N_ =M, 
and Theorem l2.5l implies that the velocities of the M incoming solitons are equal to those of 
the M outgoing solitons. However, we show in the next section that these (resonant) M- 
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soliton solutions are different from the ordinary (nonresonant) multi-soliton solutions of the 
KP equation. 

We remark that Theorem l2 . 5 I determine s the complete structure of asymptotic patterns of 
the solutions u(x,y,t) given by Jl .31 for the Toda lattice equation. In the case of the ordinary 
multi-soliton solution of the KP equation, the tau-function dl.4> does not contain all the 
possible combinations of phases, and therefore the theorem should be modified. However, 
the key idea for the asymptotic analysis of using the levels of intersection is still applicable. 
In fact, one can find from the same argument that the asymptotic velocities for the ordinary 
M-solitons are given by C2i-i,2i = fe-i + where the TM-function is the Wronskian il.4l 
with fi = e 62 '-' + e 02 ' for i = l,...,M and k\ < ki < ■ ■ ■ < ki M - Note that the velocities are 
different from those of the resonant M-soliton solution for the Toda case. 

Finally, it should be noted that the asymptotic values w\ j := —d x logXj as \ — > ±°° show 
the sorting property of the Toda lattice equation; that is, for j = l,...,N, 

d Xj ( -kj as £->°o, 

1 dx Xj-\ ' J ' J { -Kn-j+1 as ^->-oo. 

Also, one can easily show that aj — > as \t,\ — > °o, which implies the sorting behavior, i.e., 

diag(-^i,~^ 2 ,---,-^) as 
diag(-^ Ar ,-fe iV _i,...,-^ 1 ) as !=->--°o. 

Recall here that the set {X = —kj \ i = 1, . . . , A^} contains the eigenvalues of the Lax matrix L, 
with A-i > • • • > Xn as mentioned in Jl . 151 . 



3. Intermediate patterns of soliton interactions 

In this section we describe the intermediate patterns of the resonant solitons in the x-y plane. 
The key idea is to consider the pattern as a collection of fundamental resonances. The 
fundamental resonance consists of three parameters: {k\ , £2 , &3 }, that is, the case of = 3 with 
|A^_ — N+ \ = 1. Without loss of generality, let us take Af_ = 1 and = 2, i.e., a (l,2)-soliton. 
(The case of a (2,l)-soliton is obtained from the symmetry (x,y,t) — > (— x, —y, —t) of the KP 
equation, i.e., from the duality of the determinants, Ti and T2 for = 3.) Then, with k\ < ki < 
£3, the pattern of the fundamental resonance is a Y-shape graph as shown in figure Here 
and in the following we denote with [i, j] the asymptotic line soliton with c = Cjj = ki + kj. 
Notice that I(N- - 1) =/(0) = {(ih.ife)} and I(N+ - 1) = 7(1) = {(l\i,l\z) t (r^,r\ 3 )}. 

One should note that at the vertex of the Y-shape graph each index appears exactly twice 
as the result of resonance, and in figure^ those vertices form a triangle, which we refer to 
as a "resonant triangle". The resonant triangle is equivalent to the resonance condition for 
the wavenumber vectors in J2.21 Since the vertex of the Y-shape graph consists of three line 
solitons, 0, = Qj, 1 < i < j < 3, the location of the vertex is obtained from the solution of the 
equations 81 = 82 = 63, i.e., 

h-k 2 -(k\-k\) \ ( x\_(Q\-Ql-(k]-k\)t\ 

Note here that the coefficient matrix is nonsingular for k\ < ki < k 3 , and the location (x,y) 
is uniquely determined by a function of t . This implies that there always exists a Y-shape 
graph if there are three line solitons satisfying the resonance conditions J2.2b . Since the Tat + - 
function d2.3t contains all possible combinations of Af+ phases, all the vertices in the graph 
form Y-shape intersections as a result of dominant balance of three exponential terms in the 
Tjy, at each vertex. One should also note that a vertex with 4 or more line solitons is not 



On a family of solutions of the KP equation satisfying the Toda lattice hierarchy 12 




Figure 4. The Y-shape graph (left) illustrating a fundamental resonance with (ki , ki , £3 ) = 
(—1,-1,1) aIM1 me corresponding functions T)j(c), = 1,2,3 (right). The graph to the left 
represents contour lines of u(x,y,t). The circle at the level set 7(0) corresponds to the incoming 
soliton, and the stars at 1(1) correspond to the outgoing solitons. 

generic: A vertex with m distinct line solitons is obtained from the system of m equations, 
{Qi t = Qj k \ik z / z jki k=l,..., m}, in which at least m — 1 equations are linearly independent. 
Then for m > 4, this system in (x,y) is overdetermined, so that the solution exists only for 
specific choices of 0? for fixed values of t . In the cases of both ordinary and resonant 2- 
soliton solutions, the two pairs of solitons as y — * ±°° are the same, and therefore there are 
only two independent equations. Also, as mentioned before, the ordinary 2-soliton solution 
needs a balance of four exponential terms to realize an X-shape vertex. However, this balance 
cannot be dominant over a balance of three terms with the t# + -function given by d2.3l >. In 
what follows, we show that the X-shape vertex of an ordinary 2-soliton solution is blown up 
into a hole with four Y-shape vertices for the resonant 2-soliton solution. 

We now consider the case with N- = 2 and = 2, which describes the resonant 2- 
soliton solution. We can start with the graph in figure |4] having ki < Ic2 < £3. Then we 
add IC4 with £3 < £4. From Theorem l2.5l we find that both asymptotic solutions for y — » ±°° 
consist of the solitons with [1,3] and [2,4]. With k\ < Ic2 < £3 < £4, the velocity C2,4 of the 
additional soliton [2,4] as y — > — °° satisfies C2,4 > c 23 > ci,2- F° r sufficiently large negative 
values of f, the [2,4] soliton starts in the left side of the [1,3] soliton and first intersects with 
the [1,2] soliton; then the resonance condition determines that the [1,2] and [2,4] solitons 
merge and make a new outgoing soliton [1,4]. Since the N+ solitons consist of [1,3] and [2,4], 
this [1,4] soliton first branches to [1,3] and [3,4]. Then the intermediate [3,4] soliton now 
intersects with the [2,3] soliton to form the [2,4] outgoing soliton. (Note that C3.4 is the largest 
velocity among these solitons.) The process forming a resonant 2-soliton is shown in figure[5] 
Note here that there are four vertices in the interaction pattern, which correspond to the four 
resonant triangles in the c-r| plane. 

One should also note that the [2,4] soliton cannot intersect directly the [1,3] soliton unless 
a [1,2] soliton or a [3,4] soliton are created as intermediate solitons. The graph of this latter 
case is obtained from figure[5]by letting (x,y,t) — > (— x, —y, —t). Also note that the ordinary 
2-soliton solution with those same parameters (k\,... for {fi,fi} in il.6\ has different 
asymptotic solitons, namely [1,2] and [3,4], and, because of the missing exponential terms in 
the T2-function, this ordinary 2-soliton solution cannot have resonant interactions; that is, no 
resonant triangle can be formed with only those exponential terms. This is also true for any 
ordinary multi-soliton solutions of the KP equation. 
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Figure 5. A resonant 2-soliton solution u[x,y,t) (left) with (k\ .k2.k3.k4) = ( — 1, — |, |,2) and 
the corresponding functions T|/(c) (right). Both incoming and outgoing solitons correspond to 
the interstections marked by circles at the level set /(l). 



We can continue the process of adding n new incoming solitons to the graph in figure [5] 
to get a (2 + n,2)-soliton solution. One can also add m new outgoing solitons to the new graph 
to obtain a (2 + n, 2 + ra)-soliton solution. This last step can be done by adding m incoming 
solutions to a (2,2 + n)-soliton solution, which is simply obtained by the Jt rotation (i.e., 
(x,y) — > (—x, —y)) of the graph of the (2 + n, 2)-solution using the duality of the determinant. 
Then one can show the following: 

Proposition 3.1. In the generic situation, the number of holes (bounded regions) in the graph 
of the (N-,N+)-soliton solution is (N- — l)(N + - 1). 

Proof. We use mathematical induction. The case N+ = 1 corresponds to the Burgers 
equation, and it is immediate to show that the graph of the (W_,l)-soliton solution has a 
tree shape; that is, no holes (see also Ref. 1 10|). Now suppose that the (A^_,A^+)-soliton has 
(N- — 1)(^V+ — 1) holes. Add a new phase 0jv+i, with k^ + i satisfying k\ < ■ ■ ■ < < fcjy+i, 
which produces a new, fastest, incoming [N- + 1 , N + 1 ] soliton, and assume that this solution 
intersects with the [l,N- + 1] soliton, which is the slowest outgoing soliton. Then the resonant 
process of those solitons generates a [1, N+ 1] -soliton as a (2, 1 ) process, which then intersects 
with the new slowest [l,jV+2] soliton to generate an intermediate [N-+2.N+ 1] soliton. This 
intermediate soliton interacts with the second slowest outgoing soliton, the [2,A^ + 2] soliton, 
to generate [2,A^ + 3] and [A^ + 3,A^+ 1] solitons, and so on. This process is illustrated in 
figure [6] From this figure, it is obvious that there are — 1 newly created holes; that is, if 
(N-,N+) — > (A^_ + 1,N+), the number of holes increases as 

(AL - l)(N+ - 1) -> (N- - l)(N+ - 1) + (N+ - 1)=N-(N+ - 1) . 

The case of the (N-,N+ + 1) solution can be analyzed in the same way using the duality of 
the determinants. This completes the proof. □ 

One can also show the following from Proposition l3.ll 

Corollary 3.2. In the generic situation for N- +N+ = N > 3, the total numbers of intersection 
points and intermediate solitons in a (N-,N+)-soliton solution are respectively given by 
2N-N+ - N and 3N-N+ - 2N. 
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N incoming solitons 



Figure 6. A schematic diagram illustrating the creation of new holes in the resonant interaction 
process for a (JV_ + l,iV + )-soliton solution with ,/V_ +N + =N. The new soliton [JV_ +1,JV+1] 
is assumed to have a resonant interaction with the previous outgoing soliton [1,JV_ + 1]. 

Proof. By applying mathematical induction on figure^ one can easily find that the number of 
new vertices (intersection points) is 2N + — 1 and that of new intermediate solitons is 3N + — 2. 
This yields the desired results. □ 

One should compare these numbers with the case of ordinary M-soliton solution, where 
the total numbers of holes and intersection points are j(M — 1)(M — 2) and jM(M — 1), 
respectively. The resonant process blows up each vertex in an ordinary M-soliton solution to 
create a hole, so that the total number of holes in a resonant M-soliton solution is given by 

±(M-l)(M-2) + ±M(M-l) = (M-l) 2 . 

Note also that the total number of vertices in a resonant M-soliton is four times of the vertices 
of an ordinary M-soliton, i.e. each vertex is blown up to make 4 vertices with one hole. 

FigureQshows a few snapshots illustrating the temporal evolution of a resonant 3-soliton 
solution with (ki,...,k^) — (— |, — §, — i, A, I, |). This resonant 3-soliton is similar to the 
"spider-web-like" soliton solution found for the cKP equation (cf. figure 10 in Ref. |9|), even 
though the underlying equation is different in those two cases. As described in this paper, the 
behavior is determined by the structure of the tau-function which is just the sum of exponential 
functions. The tau-functions of the KP and cKP equations have the same structure for those 
solutions. 

Figure [8] shows the temporal evolution of a (4,3)-soliton solution with (k\,...,kj) = 
(—3,-2,-1,0,1,2,3). In both figure and figure [8] it can be observed that different 
intermediate solitons mediate the interaction process at different times. Also note that, for 
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Figure 7. Snapshots illustrating the temporal evolution of a resonant 3-soliton solution u(x,y,t) 
with £ 6 ) = (-§ ,-§,-5,3,3,5) and G° = ■■■ = 9° = 0: (a) t = -10, (b) t = 0, 

(c) t = 10, (d) / = 20. Note the symmetry (x,y,t) <-> (-jc, -v, —t) in (a) and (c). 

some finite values of t, the number of holes in the solution changes. However, Proposition l3.ll 
applies in the generic situation, and the total number of holes remains (N- — 1){N+ — 1), 
namely 4 holes in figure[7]and 6 holes in figure[8] In both figures, we have set all 8? = 0, so 
that all line solitons merge initially at the origin. It should be noted that even though several 
solitons might merge at the same point for some finite values of f, generically the resonant 
interactions are always among three solitons, i.e., fundamental resonances, as explained in 
this paper. 

Finally, we would like to point out that the KP equation has a large variety of multi- 
soliton-type solutions. Among those solutions, we found that, since the Xn + -function of the 
resonant (A^_,A^+)-soliton for the Toda lattice hierarchy contains all possible combinations 
of phase terms {8, | i = l,...,N}, the interaction process for these solutions results in a 
fully resonant situation. On the other hand, the ordinary M-soliton solutions display a 
nonresonant case; that is, resonant triangles representing either (2,1)- or (l,2)-solitons cannot 
be formed because of the missing exponential terms in the tau-function. One can then find 
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a partially resonant case consisting of ordinary multi-soliton interaction with the addition of 
some resonant interactions; one such example is the case having f\ = e 9 ' + e 02 + e 03 and 
fl = e 93 + e 04 for the ^-function il .41 where the ordinary 2-soliton interaction coexists with 
resonant interactions. We will report the details of the general patterns for multi-soliton-like 
solutions for the KP equation in a future communication. 
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Appendix: The Wronskian solutions of the KP hierarchy 

In this Appendix we briefly explain how the Wronskian solution ( II. 41 is obtained from the 
Sato theory (see Ref. 1 14 1 for more details). The Sato theory is formulated on the basis of a 
pseudo-differential operator, 

L = d + u 2 d~ l +m 3 3~ 2 H , 

where 8 is a derivation satisfying = 3 _1 d = 1 and the generalized Leibnitz rule, 

3 V £ (l)U?- k 8, for VEZ. 



k=0 



k J dx k 



(Note that the series terminates if and only if v is a positive integer.) Then the KP hierarchy 
can be written in the Lax form 

^- = [<B n ,L], with a„:=(x")> , (A.l) 

where (l")>q represents the polynomial (differential) part of L n in d. Here the solution of 
the KP equation II. 1> is given by u = 2u2 with t\=x, t 2 —y and f 3 = t. 
Now writing L in the dressing form, 

L = WdW~ l , with W = 1 + wid~ l + w 2 d~ 2 H , 

the KP hierarchy becomes 

— — =-B n W -tVd". (A.2) 
dt„ 

Using (IA.U . the variables m,'s can be expressed in terms of w/s; for example, 

u 2 = -wij, 

M 3 = -W2 r t + WlWl tI , 

and so on. (Here and in the following, subscripts x and t„ denote partial differentiation.) The 
equations for wj are, for example, 

Wlfy = -?-WlWls-\-Wl >xx + 2w2j., 
W 2 .t 2 = -2W2Wis+W2,xx + 2W3s, 

and so on. Here one can easily show that a finite truncation of <W , given by 

W M := 1 +wi3~' H ^WMd~ M , 



is invariant under the equation (IA.2> . For example, the W -equation with M = 1 truncation, 
i.e. W\ = 1 + w\d~ l , is just the Burgers equation, 

wify = -2wiwi vV + wi,„. (A.3) 

For the M-truncation, consider the ordinary differential equation for a function /, 

W M d M f = fM + Wl f("-V + ... + WM f = . (A.4) 

Let {fj | j = 1, . . . ,M} be a fundamental set of solutions of (IA.4> . Then the coefficient function 
wi is expressed in terms of the Wronskian for the set of those solutions, i.e., 

w\ = - — logx M , with z M = Wr(/i,...,/ M ). 

ox 

which leads to a solution of the KP equation, 

a a 2 

u = 2u 2 = -2— w\ =2 T -rlogTM- 

ox OX L 
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Recall that for M — 1 this equation gives the well-known Cole-Hopf transformation between 
the Burgers equation for w\ and the linear diffusion equation for Ti = /. One can also show 
from jA.2l that / satisfies the linear partial differential equations, 

d f d "f f 1 n 

dF n = d^> for ra = 1 ' 2 ' - 

Thus the equations for (w\, . . . ,wm) on the M-truncation are linearizable, and the behavior of 
the solutions is expected to be similar to the case of the Burgers equation. (The M-truncated 
equation is a multi-component extension of the Burgers equation |6|.) This is one of the main 
motivations of the present study. 
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